Dynamics of Crossover from a Chaotic to a Power Law State in Jerky Flow 
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We study the dynamics of an intriguing crossover from a chaotic to a power law state as a function 
of strain rate within the context of a recently introduced model which reproduces the crossover. 
While the chaotic regime has a small set of positive Lyapunov exponents, interestingly, the scaling 
regime has a power law distribution of null exponents which also exhibits a power law. The slow 
manifold analysis of the model shows that while a large proportion of dislocations are pinned in the 
chaotic regime, most of them are pushed to the threshold of unpinning in the scaling regime, thus 
providing insight into the mechanism of crossover. 

PACS numbers: 05.65.+b, 05.45.Ac, 62.20.Fe, 05.90.+m 



The Portevin Le-Chatelier effect jl|, or the jerky flow 
is a rare example where the complex spatio-temporal dy- 
namics results from the collective behavior of participat- 
ing defects, namely, dislocations 0]. When samples (usu- 
ally thin strips) of dilute alloys are subjected constant 
strain rate deformation in a window of applied strain 
rates (e a ) and temperatures, one finds repeated stress 
drops (or yield drops). Each stress drop is associated 
with the formation of a band of dislocations. At low 
strain rates the bands are static nucleating randomly in 
space. At intermediate e a , bands nucleate one ahead of 
the other in a hopping manner. At high strain rates, 
bands propagate continuously. The bands found in differ- 
ent regimes of strain rates are considered to be different 
correlated states of dislocations within the bands Q. 

A classical explanation of the jerky flow is through dy- 
namic strain aging Q • At small velocities of dislocations 
(or e a ), solute atoms have sufficient time to diffuse to 
dislocations, and pin them. Longer they are arrested, 
larger is the stress required to unpin them. As the stress 
increases, dislocations get unpinned and move fast till 
they are again pinned due to diffusing solutes and other 
pinning centers. The process repeats itself. Further, the 
competition between the time scales associated with dif- 
fusion and dislocation mobility translates, at the macro- 
scopic level, to a negative strain rate sensitivity of the 
flow stress. This intermittent sequence of loading and 
unloading, and the negative flow rate sensitivity are typ- 
ical features of many stick-slip situations such as fault 
dynamics |4, frictional sliding [5| and peeling of an ad- 
hesive tape p and charge density waves pjBI- ( Indeed, 
the similarity of the PLC effect with charge density waves 
has been studied in detain in Ref. The power law 

statistics of stress drops 0, || at high e a is similar to 
those in earthquakes Q and many other power law sys- 
tems 0, |nj which however are seen at low drives. 

This rich spatio-temporal dynamics has defied a proper 
understanding due to lack of techniques to describe the 
collective behavior of dislocations . Recent studies us- 
ing nonlinear dynamical methods have shown that a rich 
body of dynamical correlations is hidden in the stress- 
strain curves of jerky flow [l2j . More recently, an in- 
triguing crossover from a chaotic to power law regime 



has been observed as a function of strain rate 0,0. As 
the crossover is seen in single and poly-crystals, it ap- 
pears to be insensitive to the microstructure. 

This crossover is unusual in a number of ways. First, it 
is a rare example of a transition between two dynamically 
distinct states. Chaos is characterized by self similarity of 
the attractor and sensitivity to initial conditions arising 
from a few degrees of freedom. In contrast, power laws 
are scale free and are infinite dimensional. Thus, most 
systems exhibit either of these regimes. To the best of our 
knowledge, this is one of the two experimental situations 
where both are seen in the same system, the other being 
in hydrodynamic turbulence 01 ■ Second, the power law 
in the PLC effect ( as also in turbulence) is seen at high 
drive rates. On the other hand, power laws observed in 
varied systems, usually seen in slowly driven dissipative 
systems, are conventionally explained by invo king the 
concept of self-organized criticality (SOC) [ToL 111! Tl4|. 
Thus, one suspects that the dynamical features of the 
power law here to be closer to turbulence than the con- 
ventional SOC systems as we will show. 

Recently, we have succeeded in explaining this 
crossover in the dynamics of the PLC effect by extending 
the dynamical model for the PLC effect due to Anan- 
thakrishna and coworkers 0, 0] . The extended model 
also explains the three types of dislocation bands seen 
with increasing strain rate |17|. However, the dynamics 
of the crossover has not been elucidated. A natural tool 
for characterizing the crossover is to follow the Lyapunov 
spectrum as a function of strain rate. This will be sup- 
plemented by the slow- manifold approach 0, 0] which 
allows us get a geometrical picture of the changes occur- 
ring in the configuration of dislocations during crossover. 

We shall briefly describe the extended dynamical 
model. The original dynamical model |l5j captures the 
well separated time scales implied in dynamic strain ag- 
ing by using the fast mobile p m {x, t), immobile pi m (x, t), 
and Cottrell type of dislocations p c (x, t). Then, all quali- 
tative features of the effect emerge due to nonlinear inter- 
action of these populations, assumed to represent collec- 
tive degrees of freedom. In spite of the idealized nature 
of the model, it has been successful in explaining several 
generic features of the effect, notably - the occurrence of 
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the effect in a window of strain rates ( or temperatures) , 
and the emergence of negative strain rate sensitivity of 
the flow stres s Il5l ITaj . The model also predicts chaotic 
stress drops |20j which has been subsequently verified 
by analyzing experimental signals [9|, Il2| , and hence has 
the right framework to model the crossover. The equa- 
tion of motion of the scaled dislocation densities fl5lll8| 
p m (x,t), Pi m {x,t) and p c (x,t) are given by : 

pm = ~bop m — p m pim + pim — apm + <f>™ffpm 

+ (D/pim)((t>?ff(x)pm)xx, (1) 

Pirn. = b (b Pm - PmPim ~ Pirn + Clp c ), (2) 

Pc = c{p m - p c ). (3) 

The overdot and the subscript x refer respectively to the 
time and space derivatives. The first term in Eq.(l), 
6oPm, refers to the immobilization of two mobile disloca- 
tions due to formation of locks, p m pim to the annihilation 
of a mobile dislocation with an immobile one, pi m to the 
remobilization of the immobile dislocation due to stress 
or thermal activation, and ap m represents the immobi- 
lization of mobile dislocation due to the aggregation of 
solute atoms. Once a mobile dislocation starts acquiring 
solute atoms, we regard it as the Cottrell type p c . They 
eventually become immobile as more solute atoms aggre- 
gate. The aggregation of solute atoms can be regarded as 
the definition of p c , ie., p c = J_ dt' p m (t')exp[—c(t — t')]. 
(See Ref. 0.) 4>™tfp m refers to the rate of production 
of dislocations due to cross glide. This depends on the 
velocity of mobile dislocations taken to be V m (</>) = 4>™f f , 
where <f) e ff = (0— h^jp im ) is the scaled effective stress, <j> 
the scaled stress, m the velocity exponent and h a work 
hardening parameter. In the original model cross-slip has 
been used as a source of dislocation multiplication which, 
however, is intrinsically nonlocal. During cross-slip dislo- 
cations leave the slip plane due to, for instance, the effect 
of repulsive internal stresses and then slip back onto the 
slip plane. This mechanism transports dislocation den- 
sities through space. This is known to lead to a 'diffu- 
sive' type coupling [2|. Let Ax be an elementary length. 
Then, the flux flowing from x ± Ax and out of x 

is given by $(x) + | [ 4>(x + Ax) - 2$(s) + $(x - Ax)}, 
where $(x) = p m {x)cj)™ff{x). Here p is the probabil- 
ity of spreading into neighboring elements. Expanding 
&(x ± Ax) up to the leading terms, we get p m §™^ + 

| d a^a"^ (Ax) 2 . Since cross-slip spreads only into re- 
gions of minimum back stress arising from p im ahead of 
it, we use Aa; 2 =< Ax 2 >= f 2 p~^. Here < . > refers to 
the ensemble average and f is an elementary (dimension- 
less) length with V = pf 2 /2. The scaled constants, a, bo 
and c refer, respectively, to the concentration of solute 
atoms slowing down the mobile dislocations, the reacti- 
vation of immobile dislocations, and the diffusion rate of 
solute atoms. The orders of magnitudes of these con- 
stants are known from the basic mechanisms and their 
correspondence with experimental quantities |18| . Defin- 
ing e, d and I as the scaled strain rate, effective modulus 



of the machine and the sample, and length of the sample, 
respectively, the machine equation reads 

<P = d[e-(l/l) [ p m {x,t).<l>7fMt)dx], (4) 
Jo 

Global coupling in Eq. (4) is similar to studies on space 
charge currents in semiconductors where the integrated 
electric field balances the applied voltage. 

In the domain of the crossover parameter, e, our earlier 
analysis EHH shows that the effect is observed between 
a forward Hopf bifurcation at low strain rate and a re- 
verse one at high e. The reverse Hopf bifurcation implies 
decreasing amplitudes of stress drops as in experiments. 
All the interesting dynamics, including chaos, is seen in 
this regime as shown in Ref. 0, 0] . 

We discretize the specimen length into N equal parts, 
and solve for p m (j,t), pi m (j,t), p c (j,t), j = 1,...N, and 
4>(t). Due to the widely differing time scales, appropri- 
ate care is taken in the numerical solutions by using a 
variable step fourth order Runge-Kutta scheme with an 
accuracy of 10 -6 for all the variables. The spatial deriva- 
tive is approximated by its central difference. The initial 
values are taken as the steady state values for the vari- 
ables (as the long term evolution does not depend on the 
initial values) with a Gaussian spread along the length 
of the sample. In experiments, the ends of the sample 
have large strains induced due to high stress concentra- 
tions at the grips. To mimic the strain, we set pi m (j,t), 
j = 1 and N to values two orders higher than the rest 
of the sample. Further, as bands cannot propagate into 
the grips, we use p m (j, t) = p c (j, t) = at j = 1 and N. 
For the numerical work, we use a = 0.8, bo = 0.0005, c = 
0.08, d = 0.00006, m = 3.0, h = and V = 0.5. How- 
ever, the results hold true for a wide range of param- 
eters values in the instability domain including that of 
T>. For these values, the PLC effect is seen in the range 
10 < e < 2000. Chaotic stress drops are seen at low 
strain rates and power law statistics at high e [l6|. 

We identify the chaotic regime by calculating the Lya- 
punov exponents, Xi (i = 1,..,M = 3./V + 1), using 
Eqs. (1-4). (The various systems sizes studied from 
N = 100 — 3333 show a rapid convergence of the re- 
sults even around 300.) The largest Lyapunov exponent, 
Ai, is obtained by averaging over 15000 time steps after 
stabilization. Ai becomes positive at e = 35, reaching a 
maximum at e = 120, and practically vanishes around 
250. (See Fig. 3a of Ref. [16|. Periodic states are seen 
prior to chaos.) In the chaotic region, the distribution of 
Lyapunov exponents is quite broad. A plot for e = 120 
is shown Fig. la. Of these, only 6.2% of M(= 1051) 
are positive. As e increases to 280, concomitant with 
the decrease in the maximum Lyapunov exponent to a 
small value, ~ 5.2 x 10~ 4 , the number of null exponents 
(almost vanishing) increases gradually reaching a value 
w 0.38M in the range [-0.00052, 0.00052] (compared to 
only a few for e = 120). For e > 250, below a resolution 
~ 10~ 4 , even as the first few exponents are distinguish- 
able, most cross each other as a function of time, but the 
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FIG. 1: (a) Distribution of the Lyapunov exponents at e = 
120 for M = 1051. (b) and (c) Distribution of the null expo- 
nents at e = 280 for M = 10000. In (c) (+) refers to positive 
and (•) to negative null Lyapunov exponents. 
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FIG. 2: (a) Bent slow manifold Si and Si (thick lines) with 
a simple trajectory for e = 200 and m — 3. Inset: p m (dotted 
curve) and <j>. (b) Same trajectory in the (cj>,Pi m ,p m ) space. 



(time averaged) distribution remains unaffected. The fi- 
nite density of null exponents has a peaked nature in the 
interval 250 < e < 700. A plot is shown in Fig. lc for 
e = 280 which can be fitted to a power law -D(|A|) ~ |A| ' 
as shown in Fig. lc with 7 = 0.6. It may be pertinent 
to state here that the uncorrelated bands and hopping 
bands are seen in the chaotic regime while the continu- 
ously propagating bands are seen at high strain rate 01 • 
As yield drops are caused by unpinning of disloca- 
tions from pinned configuration, we first need to iden- 
tify these configurations. This can be done by using 
the slow manifold approach. Here we recall some rele- 
vant results [Lsl ITflj on the slow manifold of the origi- 
nal model (D = 0) for further use when D 7^ 0. Slow 
manifold expresses the fast variable in terms of the slow 
variables, conventionally done by setting the derivative 
of the fast variable to zero |18l Il9| . Here, p m = gives 
Pm = Pm{Pimi 4>)- Instead, we use p m in terms of a sin- 
gle slow variable 8 = 4> m — p lm — a. We note that 8 
takes on small positive and negative values as both pi m 
and <f> are small and positive. Using p rn = g(p m ,4>) = 
—boPm + Pm8 + Pim = 0, and noting that p m > 0, we 
get two solutions p m = [S + (S 2 + 46 /Oim) 1 ^ 2 ]/2foo, one 
for 8 < and another 8 > 0. For regions of 8 < 0, as 
bo is small ~ 10 -4 , we get p m /pi, n —1/8 which takes 
on small values defining a part of the slow manifold, 
Since physically pinned configuration of dislocations im- 
plies small mobile density and large immobile density, we 
refer to the region of S2 as the 'pinned state of disloca- 
tions'. Further, larger negative values of 8 correspond to 
strongly pinned configurations, as they refer to smaller 
ratio of Pm/Pim- Corresponding to 8 > 0, another con- 
nected piece 5*i is defined by large values of p m , given by 
p m pa 8/bo, which we refer to as the 'unpinned state of 
dislocations'. S2 and Si are separated by 8 — 0, which 
we refer to as the fold line [18J, [l9j (see below). A plot of 
the slow manifold in the 8 — p m plane is shown in Fig. 2a 
along with a simple monoperiodic trajectory describing 
the changes in the densities during one loading-unloading 
cycle. The inset shows p m {t) and 4>{t). For complete- 
ness, the corresponding plot of the slow manifold in the 
(Pm, Pim, 4>) space is shown in Fig. 2b, along with the 
trajectory and the symbols. Note that S2 and S\ are 
separated by 8 = <f> m — pi m — a = 0, and hence the name 



fold line. In Fig. 2a, as the trajectory enters S2 at A and 
moves into S2, 8 takes a maximum negative value at B. 
Then 8 increases as the trajectory returns to A' before 
leaving £2. The corresponding segment is ABA' in Fig. 
2b, which is identified with the flat region of p m (t) in the 
inset of Fig. 2a. As the trajectory crosses 8 = 0, dg/dp m 
becomes positive and it accelerates into the shaded region 
(Fig. 2a) rapidly till it reaches p m = 8/2bo. Thereafter it 
settles down quickly on S± decreasing rapidly till it reen- 
ters S2 again at A. The burst in p m (inset in Fig. 2a) 
corresponds to the segment A' DA in Fig. 2a and b. 

Now consider the stress changes as the state of the sys- 
tem goes though a burst of plastic activity. For T> = 0, 
Eq.(4) reduces to <p = d[e — e p ], where e p = <f) m p m de- 
fines the plastic strain rate. Since p m is small and nearly 
constant on S2, stress increases monotonically. However, 
during the burst in p m [A' DA in the inset), e p (t) exceeds 
e leading to an yield drop. Since p m grows outside S2, 
8 = separates the pinned state from the unpinned state, 
and hence 8 = physically corresponds to the value of 
the effective stress at which dislocations are unpinned. 

Now, we extend the slow manifold analysis to the case 
when spatial coupling D ^ to study the changes in 
the spatial configuration of dislocations as we go from 
the chaotic to scaling regime. Using h = 0,<p e ff = 
<t>, the plastic strain rate e p (t) is given by e p (t) = 

< l >m i t )j I Pm(x,t)dx = c/> m (t)p m (t), where p m {t) is the 
mean mobile density (= J2j Pm(j, Thus, the yield 

drop is controlled by the spatial average p m (t) and not 
by individual values of p m (j)- Since the yield drop oc- 
curs when p m (t) grows rapidly, it is adequate to examine 
the spatial configurations on the slow manifold at the on- 
set and at the end of typical yield drops. Figures 3 a,b 
and 3 c,d show respectively, plots of j, 8(j), p m (j) for the 
chaotic state e = 120 and the power law state e = 280, at 
the onset and at the end of an yield drop. It is clear that 
for e = 120, both at the onset and at the end of a typical 
large yield drop (Fig. 3a, b), most p m (j)s are small with 
large negative values of S(j), i.e., most dislocations are 
in a strongly pinned state. The arrows show the increase 
in Pmii) at the end of the yield drop. In contrast, in the 
scaling regime, for e = 280, most dislocations are at the 
threshold of unpinning with S(j) ~ 0, both at the onset 
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FIG. 3: Dislocation configurations on the slow manifold at 
the inset and at the end of an yield drop: (a) and (b) for 
e = 120 (chaos), and (c) and (d) for e — 280 (scaling). 

and end of the yield drop (Fig. 3c, d). This also implies 
that they remain close to this threshold all the time (Fig. 
3d). Since 5(j) w (for most j's) refers to a marginally 
stable state, it can produce almost any response. This in 
turn implies that the magnitudes of yield drops A(f> are 
scale free. We have verified that the edge-of-unpinning 
picture is valid in the entire scaling regime for a range of 
N = 100— 1000. Further, the number of spatial elements 
reaching the threshold of unpinning 6 = 0, during an 
yield drop increases as we approach the scaling regime. 

Several comments may be in order on the dynamics of 
the crossover. First, the crossover itself is smooth as the 
changes in the Lyapunov spectrum are gradual. Second, 
the power law here is of purely dynamical origin. This 
is a result of the reverse Hopf bifurcation at high strain 
rates which limits the average stress drop amplitude to 
small values 0|. Third, our analysis shows that 
the power law regime of stress drops occurring at high 
strain rates belongs to a different universality class as it 
is characterized by a dense set of null exponents. As zero 



exponents correspond to a marginal situation, their finite 
density physically implies that most spatial elements are 
close to criticality. This is supported by the geometrical 
picture of the slow manifold where most dislocations are 
at the threshold of unpinning, 5 = 0. In contrast, going 
by the few reports, the marginal nature of conventional 
SOC models does not display any characteristic feature in 
the Lyapunov spectrum 21]. (For instance, no zero and 
positive exponents, a positive and zero exponent, zero 
exponent in the large N limit etc have been reported 
|21j.) More significantly, the dense set of null Lyapunov 
exponents themselves follow a power law. Further, we 
note that the Lyapunov spectrum evolves from a set of 
both positive and negative, but few null exponents in the 
chaotic region, to a dense set of marginal exponents as 
we reach the power law regime. Thus, the dense set of 
null exponents in our model is actually similar to that 
obtained in shell models of turbulence where the power 
law is seen at high drive values |2^|. However, there are 
significant differences. First, we note that the shell model 
|22j ] cannot explain the crossover as it is only designed 
to explain the power law regime. Further, the maximum 
Lyapunov exponent is large for small viscosity parameter 
(Ai oc viscosity^ 1 / 2 ) in shell models (2^ in contrast to 
near zero value in our model. 

In conclusion we have demonstrated that the changes 
in the Lyapunov spectrum provides a good insight into 
the dynamical mechanism controlling the crossover. The 
slow manifold analysis, applied for the first time to study 
the crossover, is particularly useful in giving a geometri- 
cal picture of the spatial configurations in the chaotic and 
scaling regimes. This picture explains the origin of small 
amplitude stress drops at high strain rates. Finally, as 
far as we know, this is first fully dynamical model which 
exhibits a crossover from chaotic to power law regime and 
should be of interest to the area of dynamical systems. 

This work is supported by Department of Science and 
Technology, New Delhi, India. 
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